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ABSTRACT 

We have used the SHARC II camera at Caltech Submillimeter Observatory 
to make 350 /im and 450 fim images of the Vega dust disk at spatial resolutions 
(FWHM) of 9'.'7 and ll'.'l, respectively. The images show a ring- like morphology 
(radius ~ 100 AU) with inhomogeneous structure that is qualitatively different 
from that previously reported at 850 /zm and longer wavelengths. We attribute 
the 350/450 /zm emission to a grain population whose characteristic size (~ 1 
mm) is intermediate between that of the cm-sized grains responsible for emission 
longward of 850 fim and the much smaller grains 18 /im) in the extensive 
halo, visible at 70 /im, discussed by Su et al. (2005). We have combined our 
submillimeter images with Spitzer data at 70 /^m to produce 2-d maps of line- 
of-sight optical depth (relative column density). These "tau maps" suggest that 
the mm-sized grains are located preferentially in three symmetrically-located 
concentrations. If so, then this structure could be understood in terms of the 
Wyatt (2003) model in which planetesimals are trapped in the mean motion 
resonances of a Neptune-mass planet at 65 AU, provided allowance is made for the 
spatial distribution of dust grains to differ from that of the parent planetesimals. 
The peaks of the tau maps are, in fact, located near the expected positions 
corresponding to the 4:3 resonance. If this identification is confirmed by future 
observations, it would resolve an ambiguity with regard to the location of the 
planet. 
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1. Introduction 

The observable behavior of circumstellar dust particles under the influence of gravity 
and radiation pressure can provide information on the locations and masses of unseen planets 
(see, for example, Wyatt (2003, 2006); Deller & Maddison (2005)). A particularly suitable 
object for study is the AO star Vega at 7.76 pc whose disk is seen nearly face-on, as suggested 
by the low inclination (5°) of the stellar rotation axis (GuUiver, Hill & Adelman 1994) . Images 
at 850 //m and 1.3 mm (Holland et al. 1998; Wilner et al. 2002; Koerner, Sargent & Ostroff 
2001) indicate a partial dust ring (r ~ 100 AU) dominated by two unequal clumps interpreted 
as collisional debris from resonantly-trapped planetesimals. More recent observations with 
the Multiband Imaging Photometer for Spitzer (MIPS) (Su et al. 2005) have revealed an 
extensive halo representing small (~ 2 fim) and medium-sized (~ 18 fim) grains blown out 
from the ring by radiation pressure, the ring itself being dominated by large grains of radii 
a > 180 iim. 

Dynamical modeling shows the spatial distribution of the ring to be consistent with the 
2:l(u) resonance of a Neptune-mass planet which has undergone migration (Wyatt 2003). 
Such a model predicts the existence of other populated resonances (for example 3:2 and 
4:3), the detection of which would help confirm the model and provide new constraints. We 
present new observations at 350 /im and 450 /im which bear on this issue. 

2. Observations and Data Reduction 

Vega was observed at 350 fim and 450 fira with the SHARC II camera (Dowell et 
al. 2003) at the Nasmyth focus of the Caltech Submillimeter Observatory (CSO) on UT 

2004 Sep 17-19 (350 /xm; r225GHz ^ 0.041), 2005 Apr 23 (350 /xm; r225GHz ^ 0.036), and 

2005 Jun 11-13 (450 //m; T225ghz — 0.048), where T225ghz represents the zenith value of 
atmospheric optical depth. The on-source total integration times were 5.0 hr, 2.7 hr and 9.7 
hr, respectively. Similarly to our Fomalhaut observations (Marsh et al. 2005), the telescope 
was scanned in an oscillatory fashion in azimuth and elevation with peak-to-peak amplitudes 
of 30" — 100". The Vega images were generated with an iterative code ( "sharcsolve" ) similar 
to CRUSH^. To reduce spatial 1// noise in the center of the image, the sky intensity was 
forced to be zero beyond 35" from the star (analogous to spatial chopping with amplitude 
70" in all directions) , except for the final iteration which allowed non-zero intensities outside 
35" radius for noise-evaluation purposes. After the application of Gaussian smoothing (with 
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a kernel of 5'.'5 FWHM), the angular resolutions (FWHM) of the output images were ^'.7 
and ll'.'l, at 350 jim and 450 /xm, respectively. Absolute calibration was accomplished with 
hourly, interspersed observations of point sources, and is based on assumed Neptune fluxes 
of ^350 = 92.2 and 86.5 Jy in 2004 and 2005, respectively, and 5*450 = 65 Jy. We estimate 
Icr calibration accuracies of 30% and 1" in absolute flux and pointing, respectively. The flux 
error includes the effect of subtraction of the slowly- varying background (~ 5% level over 
spatial scales of a few tens of arcseconds at both wavelengths) . 

The point source response function (PSF) was obtained from observations of Mars at 
350 /xm and Neptune at 450 /um, rotated and coaddcd to recreate the rotational smearing 
due to changing parallactic angle. The nonzero angular diameters of those planets (3'.'5 and 
2'.'3, respectively) broadened the estimated PSFs by 6% and 2%, respectively; the effect on 
subsequent processing was found to be negligible within prevailing errors. 

Figure 1 shows the observed images at 350 fim. and 450 /xm, before and after subtracting 
the estimated photospheric contributions of 35 mJy and 21 mJy, respectively. The integrated 
flux densities of Vega (before photospheric subtraction), within a circular aperture of 30" 
radius, at 350 /im and 450 //m are 500 ± 150 mJy and 150 ± 45 mJy, respectively. 

The observed images show clearly the ring morphology of the Vega disk, and suggest 
inhomogeneous structure. Figure 2 shows the 350 //m intensity variation as a function of 
azimuth, calculated in an annulus of inner and outer radii 6'.'9 and 13'.'9, respectively. Three 
peaks are apparent, at azimuths of —90°, 0°, and 120°, with amplitudes of 3(7, 4(j, and 2a 
above the local mean level respectively, where cr represents the statistical measurement noise. 

3. Mapping the Relative Column Density of Dust 

We have mapped the line-of-sight optical depth (relative column density of dust) in the 
Vega disk as a function of 2-d location in the disk plane using the DISKFIT procedure (Marsh 
et al. 2005), assuming in this case a geometrically-thin face-on disk. The output, referred 
to as a "tau map," is estimated using a set of observed images at multiple wavelengths, 
taking full account of the corresponding PSFs. It is assumed that the local temperature 
of each dust grain component is determined by the energy balance of individual grains in 
the stellar radiation field using the results of Backman & Paresce (1993). The current 
version of the code makes simultaneous estimates of the tau maps corresponding to the 
different grain components, each of which is characterized by the parameters Aq and 
where Aq represents the wavelength above which the grains radiate inefficiently and has an 
approximate correspondence with the grain radius, a, and (5 is the power-law index of the 
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wavelength dependence of opacity such that oc A ^. The spatial resolution of the tau 
maps exceeds that of the raw images due to implicit deconvolution of the PSFs. 

Tau maps were made using data at three wavelengths, by combining our 350 & 450 /im 
images with a 70 //m MIPS image at the fine (5") pixel scale. We employed a two-grain 
model consisting of the 18 //m grains which dominate the 70 //m emission (Su et al. 2005) 

and a population of larger grains whose size was chosen based primarily on the observed 
spectral slope between 350 ^ra and 70 /um. Specifically, assuming 13 = 1 (Dent et al. 2000), 
the spectral slope implies Tdust ~ 50 K for the larger grains; if we further assume a size- 
temperature dependence corresponding to the grain composition (silicate-carbon mix) used 
by Su et al. (2005) in their Vega modeling, we obtain a > 100 /^m. Since this limit was based 
on the extremes of the flux error bars and of the relative contributions of ring and halo at 
70 //m, we have adopted a larger value (1 mm) as being more likely to be representative of 
the grains reponsible for 350/450 //m emission. The results are presented in Figure 3, which 
shows tau maps for the two grain components (18 /xm and 1 mm) separately and also the 
total optical depth (lower panel). Superposed on the latter, for comparison, are the positions 
of previously-reported emission peaks from longer wavelength data at 850 /xm and 1.3 mm 
(Holland et al. 1998; Koerner, Sargent & Ostroff 2001; Wilner et al. 2002). 

The tau map for 1 mm grains shows inhomogeneities, the reality of which we have 
assessed using a test based on fits to the 350 /im and 450 /xm data. Specifically, an 
azimuthally-uniform ring resulted in a 10% increase in relative to an unconstrained tau 
map; taking into account the number of independent data points (255), this translates into 
a relative probability ~ 5 x 10~^ corresponding to a 4.6(j deviation. We therefore conclude 
that the azimuthal structure we see in the ring is at the 4(7-5(7 significance level, consistent 
with our findings based on the S/N oi peaks in the observed 350 //m image (see Figure 2). 

The tau map for 18 ^m. grains is dominated by a ring of larger scale (r ~ 140 AU) 
than for the 1 mm grains (r ~ 100 AU). There is also lower-level diffuse structure masked 
by surrounding noise spikes, the latter resulting from the radial increase in tau estimation 
error corresponding to decreasing dust temperature. We can, however, smooth out these 
fiuctuations with azimuthal averaging to show the underlying diffuse structure. Figure 4 
shows the result in the form of radial profile plots of the 1 mm grain component (filled 
circles) and the 18 /im grain component (open circles). The error bars reflect the tau map 
estimation errors after a \fN reduction from azimuthal averaging in annuli of width 17.5 
AU. It is apparent that for r > 200 AU, the radial falloff in column density of 18 ^m. grains 
is consistent with the 1/r variation reported by Su et al. (2005). Inside 200 AU, however, 
some of these grains may be concentrated in the the ring structure, but the proportion is 
subject to model uncertainties related to the larger grain component. 
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4. Interpretation 

The observations show quahtatively different appearance of the Vega dust ring at 
350/450 fim and > 850 /im. Furthermore, the 850 fim flux (23 mJy) predicted by our 
tau map accounts for only a fraction of either of the pubhshed values of 850 /xm flux, which 
are 45.7 ± 5.4 Jy (Holland et al. 1998) and 91 mJy (Holland et al. (2005), as quoted by 
Su et al. (2005)). The observations can, however, be reconciled using a model in which the 
ring emission at 350/450 /im is dominated by 1 mm grains with a opacity law, and the 
emission longward of 850 /xm is dominated by an additional component of much larger grains 
with flat spectral opacity (/? ~ 0); the latter would require grain sizes of at least a centimeter 
(Pollack et al. 1994). 

The morphology of the 850 /xm image (with its single dominant peak at ~ 45° in position 
angle) has been attributed to dynamical effects, the most important of which is the 2:l(u) 
resonance of a Neptune-mass planet 65 AU from the star (Wyatt 2003). Those calculations 
also show that other resonances are present, principally 3:2 and 4:3, the latter of which shows 
a 3-fold spatial symmetry consistent with what we see at 350/450 fim. In this regard we note 
that for a planet of given mass, the requirement that the dust grains not be dislodged from 
resonance by radiation pressure sets a lower limit, Ocrit, to the grain size (Wyatt 2006). For 
a Neptune-mass planet orbiting Vega (luminosity 60Lq; mass 2.5M0), Wyatt's expression 
leads to Ocrit = 900 /im, consistent with our choice of a = 1 mm for the tau maps. 

Based on the 850 /xm data alone, an ambiguity exists for the orbital direction of the 
planet, shown counterclockwise in Figure 15(b) of Wyatt (2003). But if our identification of 
the 350/450 iim structure with the 4:3 resonance is correct, then the ambiguity is resolved; 
the planet is then located diametrically opposite to the star from the SE clump in our Figure 
3 and its orbit is clockwise. The predicted locations of the 4:3 resonance peaks are then given 
by the arrows in Figure 2, which are in substantial agreement with observation. 

It would then remain to explain why the cm-sized grains (seen at 850 /im) and mm-sized 
grains (seen at 350/450 jim) exhibit two different resonance patterns. We note, however, that 
the analysis by Su et al. (2005) suggests that the dust has been released during a relatively 
recent event, most likely the collision of two planetesimals. This being the case, our results 
could be understood if the two planetesimals involved were located in the 4:3 resonance, 
and that the collision released smaller particles than the ones which were there before, thus 
affecting the 350/450 fim appearance much more than that at 850 /xm. The 4:3 resonance is, 
in fact, a relatively likely place for a collision to occur, since the collision timescale decreases 
with decreasing orbital semimajor axis, and hence is shorter in the 4:3 resonance than in 2:1. 

We have simulated this scenario with dynamical modeling based on a Neptune-mass 
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planet at 65 AU, assuming orbital parameters similar to Wyatt (2003), and further assuming 
that the 4:3 resonance is populated by smaller grains (1 mm) than the other resonances (1 
cm). The results are presented in Figure 5. Comparison of the predicted 350 fim and 850 /xm 
images with the observations (Figure 1 of this paper and Holland et al. (1998), respectively) 
shows that the model successfully accounts for the wavelength dependence of the observed 
structure. Our analysis therefore suggests that different resonances can be populated by grain 
populations of different sizes, depending on the coUisional history of the planetesimals. This 
means that the dust distribution does not necessarily mimic the planetesimal distribution 
as assumed by Wyatt (2003). It will thus be important to confirm the detection of the 
4:3 resonance and constrain the grain sizes involved using future observations at multiple 
wavelengths. 

This work was performed by the Jet Propulsion Laboratory, Cahfornia Institute of Tech- 
nology, under contract with the National Aeronautics and Space Administration. Research 
at the Caltech Submillimeter Observatory is supported by NSF grant AST-0229008. 
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Fig. 1. — Observed and photosphere-subtracted images of Vega at 350 fim and 450 fim. 
The intensity scale (shown by the horizontal bar at the top of each image) is in units of 
mJy arcsec"^ and the orientation is such that north is up and east is to the left. The RMS 
measurement noise on the images is 0.032 and 0.026 mJy arcsec"^ at 350 fim and 450 fim, 
respectively. 
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Fig. 2. — Intensity in the Vega disk at 350 /im wavelength as a function of azimuth around 
the ring (measured east from north), both with and without subtraction of the estimated 
photosphcric contribution. Error bars represent the standard deviation of measurement 
noise. The arrows represent the expected locations of the 4:3 resonance peaks for the dy- 
namical model of Wyatt (2003) in which dust is trapped in the mean motion resonances of 
a suspected planet. 
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Fig. 3. — The estimated line-of-sight optical depth of the Vega disk at a reference wavelength 
of 350 /im, based on observed images at 350 fim, 450 fim, and 70 /im, assuming a two- 
component dust grain model in which the grain radii are 1 mm and 18 /xm. The upper 
two panels show the optical depth distribution for the two grain components separately, 
while the lower panel shows the sum of both components. Also shown on the lower plot are 
the locations of previously reported emission peaks from longer wavelength data as follows: 
(filled triangle) 850 /im— Holland et al. (1998); (x) 1.3 mm— Wilner et al. (2002); (+) 1.3 
mm — Koerner, Sargent & Ostroff (2001). 
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Fig. 4. — Radial profiles of estimated optical depth for the two grain components, charac- 
terized by grain radii of 1 mm (filled circles) and 18 /xm (open circles). 
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Fig. 5. — Results of dynamical modeling of the Vega dust disk showing the predicted intensity 
distributions of dust emission at 350 /im and 850 /im. Following Wyatt (2003) we assumed 
a Neptune-mass planet at 65 AU, and further assumed that the planet orbits clockwise, is 
currently located NW of the star and that its 4:3 resonance is populated by smaller grains 
(1 mm) than the other resonances (1 cm). Upper two panels: The 350 /im emission resulting 
from (a) the combined effects of the 2:l(u) and 3:2 resonances, and (b) the 4:3 resonance. 
The planet is indicated by an encircled cross. Lower two panels: Predicted observational 
images at (c) 350 /xm and (d) 850 /im. Images (c) and (d) were produced by adjusting 
the relative amounts of material in the two grain populations (1 mm and 1 cm) to provide 
the best match to the observed images; this required equal amounts of material (by cross- 
sectional area) in each. In this model, the 1 mm grains (4:3 resonance) contribute 80 ± 16% 
of the observed flux at 350 /im and 60 ± 12% at 850 /im, while the 1 cm grains (2:l(u) & 3:2 
resonances) contribute 20% at 350 /zm and 40% at 850 /im. All of the above images have 
been smoothed to the spatial resolution of the observations, corresponding to FWHM of 10" 
at 350 /im (CSO/SHARC II; present paper) and 14" at 850 /im (JCMT/SCUBA; Holland 
et al. (1998)). On the linear pseudocolor intensity scale, white corresponds to peak intensity 
for each individual image. 



